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For quantum theories with a classical limit (which includes the large N limits of 
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typical field theories) , we derive a hierarchy of evolution equations for equal time 
correlators which systematically incorporate corrections to the limiting classical 



evolution. Explicit expressions are given for next-to-leading order, and next-to- 



next-to-leading order time evolution. The large N limit of iV-component vector 
models, and the usual semiclassical limit of point particle quantum mechanics 
are used as concrete examples. Our formulation directly exploits the appropri- 
ate group structure which underlies the construction of suitable coherent states 
and generates the classical phase space. We discuss the growth of truncation 
error with time, and argue that truncations of the large- iV evolution equations 
are generically expected to be useful only for times short compared to a "deco- 
herence" time which scales like N 1 / 2 . 
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I. INTRODUCTION 



The time evolution of quantum systems away from equilibrium is of interest in many applica- 
tions including, but certainly not limited to, phase transition dynamics, inflationary reheating, 
and heavy ion collisions. Large JV expansions have provided a widely used technique for study- 
ing equilibrium properties in statistical physics and field theory (T|-§|, and it is natural to apply 
a similar strategy for studying non-equilibrium problems. The large JV limit (as typically for- 
mulated) is actually a special type of classical limit M . Suitable observables behave classically 
and the quantum dynamics reduces to classical dynamics on an appropriate phase space. 

Considerable work has been done examining the dynamics of far from equilibrium states 
in a variety of applications using leading large- JV time-evolution PHIoH- A major virtue of 
large iV techniques (compared to alternative wholly uncontrolled approximation schemes) is 
that one should be able to improve the approximation by systematically including sub-leading 
effects suppressed by powers of 1/JV. For a variety of equilibrium problems (such as critical 
phenomena) , this approach can work quite well |TT|-[13| . 



For initial value problems, in which one would like to choose a non-equilibrium initial state 
and then examine the subsequent time evolution, traditional formulations of large- JV expansions 
using graphical or functional integral techniques || are very awkward. A major difficulty with 
these approaches is that they generate integral equations which are non-local in time when 
sub- leading 1/JV corrections are retained. For practical (numerical) applications, one would 
vastly prefer a formulation in which locality in time is always preserved. 

In this paper, we describe a formulation of large N (or semi-classical) dynamics which leads 
to a coupled hierarchy of time-local evolution equations for equal time correlation functions. 
Our approach directly exploits the appropriate group structure underlying the construction of 
suitable coherent states and the existence of the classical limit [[| . We specifically focus on the 
time evolution of initial states chosen to equal one of these coherent states. We will give explicit 
next-to-leading order (NLO), and next-to-next-to- leading order (NNLO), expressions for the 
required evolution equations. Somewhat related hierarchies of evolution equations have been 
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discussed in several recent papers ||14|jl5| . Because of our exploitation of the underlying group 
structure, the formulation we derive is more efficient, in the sense that it requires integration 
of fewer coupled equations at a given order in 1/N. 

A major question which we discuss, but do not fully resolve, is the propagation of errors 
induced by truncating the exact (infinite) hierarchy at a given order in 1/N. It is known 
that the N — > oo limit is not uniform in time. For example, in typical large N field theories 
the characteristic time scales for scattering or thermalization are known to scale as N to 
some positive power.Q For a fixed time interval t, results obtained by integrating evolution 
equations truncated at, for example, next-to-leading order, will have only order 1/N 2 errors. 
For sufficiently large N, and fixed t, including successively higher orders in the 1/N hierarchy 
will yield more accurate results. But for fixed N and some given truncation of the 1/N hierarchy, 
it should be expected that the truncation error will grow with increasing time and eventually 
become order unity. A key question is how this "breakdown" time scales with N and the order 
of the truncation. One might hope that a next-to-leading order approximation would be useful 
[that is, have at most 0(1/N) global errors] for times of order N, while a next-to-next-to- 
leading order scheme would be useful out to times of order N 2 , etc. But it is quite conceivable 
that errors in an order-fc truncation will grow with time like (t a /N) k for some positive a, which 
would imply that all truncations break down after a time of order N 1 ^ . This behavior, which 
we consider likely, may well depend on the specific theory and choice of initial state. Available 
numerical work, such as |T^JT5|| , sheds little light on this issue. We discuss several examples 
where it is possible to argue that quantum "decoherence" produces exactly this type of limit 
on the range of validity of large N truncations. 

The paper is arranged as follows. The general framework which allows us to treat many 
theories with a classical limit in a uniform fashion is outlined in section II. This material 
is largely taken from Ref. M. Section III describes the particular class of operators we will 



1 See, for example, the end of section III of Ref. ml 
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consider, and examines the structure of their coherent state equal time correlators. Section IV 
presents the resulting time evolution equations and discusses error propagation. These general 
results are applied to the examples of point particle quantum mechanics, and a general N- 
component vector model, in section V. For point particle quantum mechanics, we argue that 
the decoherence time generically scales as Ti~ 1 ^ 2 , while for vector models it should scale as A 1 / 2 . 
A brief concluding discussion follows. 



II. COHERENCE GROUP AND COHERENT STATES 

The following slightly abstract framework is applicable to typical large iV limits (including 
O(N) or U(N) invariant vector models, matrix models, and non-Abelian gauge theories), as 
well as the h — > limit of ordinary quantum mechanics [|J. 

Consider a quantum theory depending on some parameter \ (such as % or 1/N). The Hilbert 
space (which may depend on x) will be denoted TL X . The quantum dynamics is governed by a 
Hamiltonian which we will write as (h/x) H x . This rescaling of the Hamiltonian will prove to 
be convenient, and makes the Heisenberg equations of motion take the form 

The following assumptions are a set of sufficient conditions implying that the x ~~ > limit 
is a classical limit. 

Assume there is a Lie group G (called the coherence group) which, for every value of x-, 
has a unitary representation on H x , G x = {D x {u) : u G G}. The states generated by applying 
elements of the coherence group to some (normalized) base state |0) x G 7i x , 

\u) x = D x (u)\0) x , ueG, (2) 

are called coherent states. The coherence group acts on these states in a natural way, 
D x (u') \u) x = \u'u) x . 

We assume that the coherence group G x acts irreducibly on the corresponding Hilbert 
space 7i x . In other words, no operator (except the identity) commutes with all elements of the 
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coherence group. This condition automatically implies that the set of coherent states form an 
over-complete basis for the Hilbert space TC X . It also implies that any operator acting on Ti x 
may be represented as a linear combination of elements of the coherence group. 

For any operator A acting in H x , we define its symbol A x {u) as the set of coherent state 
expectation values, A x {u) = (u\A\u) x , u G G. We assume that the only operator whose symbol 
vanishes identically is the null operator. Thus, distinct operators have different symbols, which 
means that any operator can, in principle, be completely reconstructed solely from its diagonal 
matrix elements in the coherent state basis. 

Classical observables will be associated with operators that remain non-singular as \ goes 
to zero, that is, whose coherent state matrix elements, (u\A\u') x / (u\u') x , do not blow up as 
X — > for all u, v! G G. Such operators are called classical. 

Two coherent states \u) and \v!) are termed classically equivalent (we will write u ~ u') 
if in the x ~ > limit, one can not distinguish between them using only classical operators, 
i.e., lim x ^ ^x( M ) = ^ m x^o J ^x( u ') f° r a ^ classical operators A. We assume that the overlap 
between any two classically inequivalent coherent states decreases exponentially with 1/x in 
the x ~~ * limit. 

Under these assumptions, one may show that the % limit of this theory truly is a 
classical limit §]. The assumptions hold for O(N) or U(N) invariant vector models, matrix 
models, and gauge theories Q . The quantum dynamics reduces to classical dynamics on a phase 
space T given by a coadjoint orbit of the coherence group. Formally, points in T correspond 
to equivalence classes of coherent states, r = {[u] : u G G}, with [u] = {u f G G : u ~ u'}. The 
symplectic structure on the phase space is completely determined by the Lie algebra structure 
of the coherence group. The classical Hamiltonian is just the x ~ > limit of the coherent state 
expectation of the quantum Hamiltonian, 

h c \{u) = lim H x {u) . (3) 

To have sensible classical dynamics this limit must exist, i.e., H x must be a classical operator. 
(This is why it was convenient the rescale the Hamiltonian by h/x-) The classical action is 
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S d [u{t)} = lim / dt(u(t)\i X d t -H x \u{t)}x- (4) 

Both the classical Hamiltonian (||) and the action ([|) depend only on the equivalence class of 
the coherent state |it),Q and thus do define sensible dynamics on the classical phase space. 

The preceding discussion is just a formalization of the usual picture of a classical limit. A 
quantum mechanical wave packet, with a width of order x 1//2 > behaves classically in the x ~^ 
limit, and may be associated with a point in the classical phase space. The equations of motion 
that govern the classical dynamics are just coherent state expectations of the original quantum 
evolution equations. 

III. COHERENT STATE EXPECTATIONS 

As noted earlier, the irreducibility of the coherence group implies that all operators may be 
(formally) constructed from the generators of the coherence group. Consequently, for charac- 
terizing the structure, and time evolution, of any state, one may focus attention on equal-time 
expectation values of products of coherence group generators. 

Let g denote the Lie algebra of the coherence group G. Let {ej} be a basis of g. The 
commutator of basis elements defines the structure constants, [e i: ej] = if^ek- The generators 
Ci themselves are not classical operators, but rather are 1/x times classical operators. For 
convenience, let Xi denote the rescaled generator which is a classical operator, Xi = x e i- 

Consider the coherent state expectation value of the monomial x^x^ ■ ■ - Xi k . We would like 
to find an expansion of this expectation value in powers of x- A convenient representation for 
our purposes involves subtracted expectations^ 

9iJ 2 --ik = (O^i — ^i) ' ' ' — x ik)) i (5) 

2 For the action, this is true up to temporal boundary terms which do not affect the dynamics. 

3 To simplify notation, we will omit the superscript "(fc)" when this can cause no confusion; for example, we 

(2) 

will write gij for g\ , etc. The same remark applies to the connected expectations discussed below. 
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FIG. 1. Expansion of the expectation a product of generators in terms of (a) subtracted, and (b) connected 
expectations. The shaded bubbles on the left denote full expectations of the product of generators. Each line 
with a dot on the end represents the expectation of a single generator. In (a), dashed-line bubbles correspond 
to subtracted expectations of strings of generators, while in (b) solid-line bubbles represent connected s^ k ' 
expectations. Ellipses (• ■ •) denote ordered permutations of the preceding diagram. 

where (■ • •) denotes an expectation in some coherent state, and X{ = (x^ are the expectations 
of the rescaled generators X{. Subtracted and un-subtracted expectations are related by 

? (3) 

(h,h) X ^ X ^ (h,h,h) x h x h x h x h ' ' ' x ik 

where the n-tuples (h,h, ■■■Jn) are ordered subsets of {zi, (There is no g^ 1 ' term since 

9i = (xi - = 0.) 

Alternatively, one may expand in terms of connected expectations, 

S iJ--i k = ( X H ' ' ' X ik) ■ (7) 

The difference, illustrated graphically in Fig. [I], is that expansions in terms of connected ex- 
pectations involve products of all possible 'contractions', while the terms in the expansion in 
subtracted expectations have only one string of generators 'contracted'. The difference between 
subtracted and connected expectations first arises with four generators. Explicitly, 

/ — — \ _ , (2) , (2) , (2) (2) (2) (2) 

\XiXjXkXi) — XiXjX^Xi -+- XiXj g^i -\- X^X}.g^ ~\~%iXigjk i~ x k%igij < x j x l9ik ' x j x k 9il 

+ x i gfki + x 3 9m + x k gfjl + xi gf) k + gf kl (8a) 
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(3) (3) (3) (3) (2) (2) (2) (2) (2) (2) (4) / R , x 

i X i S jkl ' X j S ikl ' X fc S ij7 + X « S ijk ' S ij S kl ' S ik S jl ' S il S jk ' S ijkl • \°°) 

The coherent state overlap (u\u') x is the generating functional for expectations of products 
of generators, since variations of the coherent state v! can bring down any desired generator 
of the Lie algebra, 5i \u') = ei\u'). The logarithm of this overlap is therefore the generating 
functional for connected expectations. By assumption, ln(M|M') x is 0(l/x) as x — > 0. This 
immediately implies that the k-th. order connected expectation is 0(x k ~ 1 )- Note also that 
the commutator of functional derivatives is the functional derivative in the direction of the 
commutator, 

/. . . t-t . . . .\conn _ / - - \conn _ / JA. A .1 . . \conn _ • fm / - . . \conn /q\ 
\ J/j^j / \ J'j^! / \ L^-d -k-jl j L A.Jij \ -^m I i \ J ) 

nr o(^) — ? ' v 1) 

By considering which connected expectations contribute to g^ k \ one may easily seeQ that 
subtracted expectations fall off roughly half as fast as the connected ones, g( 2k > ~ g( 2k ^ 1 ) = 
0(x k )- Because g( fc+2 )/g( fc ) = 0(x), expansion (^) is a power series in Xi the parameter 
measuring how close the system is to being classical. Of course, subtracted expectations may 
always be rewritten in terms of connected expectations (and vice- versa) . Ultimately, equations 
for connected^ expectations will be most useful. Nevertheless, using subtracted expectations 
as an intermediate representation is helpful because of the simple form of expansions in terms 
of subtracted expectations, as shown by Eq. @ and Eq. (|IID below. For later use, note that 

o {k) -o {k) -iYf m (r o (k ~ 2) + o {k ~ l) \ (10) 

Now consider an operator V that can (at least formally) be expanded in powers of genera- 
tors, V = Y,k Y,{i!---i k } a h-ik x h ' ' ' x ik-> f° r some se t of coefficients {a^...^}. Operators of this 

4 As gW = J2 s^ mi h^ ■ ■ • s( m ") with mi + h m„ = k and > 1 for all i, and = C(x mi_1 ), 

g(k) _ 0^mm{5J(m i -l)}^ _ ^,^min{fe-ri} y T ^ e ^ ar g est numbej- f connected diagrams n occurs when all s (mi) 

are (except for one, if k is odd, which is s*- 3 -*). 
5 Or perhaps one-particle irreducible. 
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form are well behaved for x ~~ * 0, and so are good classical operators. Using (|6]), 

f o (2) o {3) o (fe) 1 

<v) = ££ ^^1-^1+ £ — + £ + n " tfc ■ (ii) 

k {ii-ifc} [ (h,h) XllXl2 (hMM XllXl2Xl3 x h"' x h\ 

This can be packaged in an even more concise form, 

(V) = V" + ^ (W2) + & ^ lWs) + & 3h V^ 1 ^ + 0( X 3 ) (12) 

where summation on repeated indices is implied, and V = Y.kY.{i 1 ---i k } a h--i k x h " " x i k is the 
number obtained by replacing each generator in V by its coherent state expectation. Here we 
have introduced "ordered derivatives" f fa'"**) = tt— ^ — T defined by 

(fg)^- 1 ^ = f^-^g^" 1 ^ . (13b) 
i=0 

When acting on a string of generators, ordered derivatives produce a sum of products of expec- 
tations of the generators which remain after deleting the indicated generators, provided these 
appear (not necessarily contiguously) somewhere within the string in the order specified by the 
derivative. For example,[] x 2 px = 2x 2 , and xpx 2 p = x 2 p + 2xpx + px 2 = Ax 2 p. 

In the x ~~ y limit, coherent state expectations of the (rescaled) generators Xi turn into 
coordinates Xi on the classical phase space and (classical) operators acting on 7i x become 
functions on phase space, (V) = V + 0(x)- For finite Xi the successive terms in fll2|) precisely 
characterize the corrections to this classical limit. 



IV. TIME EVOLUTION 

Since operators are completely determined by their symbols, to study the time dependence 
of any observable A it is sufficient to take the coherent state expectation value of its Heisenberg 
equation of motion ([!]), 



6 If , Xi commute, then -ft — — r = \ „ — ^—L — . In this case, the ordering: does not matter, and the 

n! is needed to make up for over-counting. 



S 



> • (14) 

In other words, we assume that the initial state is precisely some coherent state \u), and wish 
to determine the subsequent time evolution. To do so, we will first find an expansion, in powers 
of Xi f° r t ne expectation of the commutator of classical operators. 



H X ,A 



A. Symbols of Commutators 

Consider classical operators A and B which (as in Section pl|) may be written as power 
series in the generators, A = J2 a h-i m ^h ' ' B = ^2 Pjvjn^h ' ' Their product is 

given by AB = J2 a h-i m Pji-jn^ii ' ' '%im.%h ' ' Using our previous result flTJ), we find 

(AB) =AB + g® (AB)V^ + g\% h (AB)^ + g\% hh (ABf^ + 0( X 3 ) (15a) 
= AB + 41 (AB)^ + sgU (AB)^ 

+ (4i sii + -a -a + -a -a) + ^ 

where now (1%, l 2 , ■ ■ •) denote ordered subsets of {ii, • • • , i m ,ji, • • • , i n }- We see from ( fL5f) that, 
to leading order, products of classical operators factorize, (AB) = (A)(B) + 0(x). 

Using the expansion (p~5|) , and the reduction formulas for operator derivatives (0), one can 
evaluate the commutator. A generic term in the result will be 



(k) 
9hh-i k 



= g^ l ^ (a^ 1 '" 1 ^ B + h A^ 11 " 2 ^ £?( /3 '"' fe ) + A^ B^ 2 '" 1 ^ -f AB^ ll '" lk ^ 

— g\ k ^ l ( B^ ll '" lk> A + ■ • • + B^ 11 ' 2 ^ A^ 3 '" 1 ^ + B^A^ 2 '" 1 ^ + BA^ ll '" lh ^ 

= ~ & l4 -i kh ) (A<MB«*<*"U - B^A^-^) 

i (J k ) _ J k ) \ ( A(hh) B {hh-h) _ B (hh) A {hh-h)\ 

+ \Phhhh-i k 9i 3 i 4 -i k hi 2 ) \si a D I\ j 

+ [& hh .., k - 9ti khhh ) (A^W"U - B^A^) + ■ ■ ■ . (16) 

The last term in the sum (|T§) is either (//z^.,., . - ///^ .„, ) (\l ^ ? -7i :/ <- '"'--j 
or (gl* j .$ j+1 -i 2j+1 ~ 9l-itX +1 h-l 3 ) (A^l^-r-hi.i) _ B (h-l 3 ) A (i,.i-h,.^ j depending on 



whether k is even or odd. Using (pi]) to reduce the differences (g(y — g^) yields the final form 
for the expectation of the commutator of classical operators. The leading term is precisely 
the Poisson bracket on the classical phase space, while subsequent terms involve successively 
higher expectations g( k >. Displaying subleading 0(x) and 0(x 2 ) terms explicitly, one finds 



-iAB]) 



} 0( X °) 



+ 
+ 



m 2 9m h + m 3 9i 2m ] (awbw - bwam 

x m (f£ h g hl4 + m 3 g hl4 + m 4 g hh )} (a^b^ - b^a^ 



m [fhh 9\ % k + fhk 9hh + fhh 9lik + fhk 9l 3 h 



+ 
+ 
+ 



j\(h)£>(hhh) _ ^{l^Aihhk) 



fhh 9mhU + fhh 9l 2 ml 4 + fhk 9l 2 h; 

f hh 9ml 2 k f lik 9l i ml 2 f l 2 l 3 9l 1 mk f l 2 k 9l 3 hm 

x ™ {fhh dhkls + fhh 3l 2 kh + fhk 9l 2 hh + fhh 9hhh 

+fl 1^2 9mlskh flih 9l 2 mkh fhk 9l 2 l s mlc, fhh 9l 2 l 3 k 

X {A^ .gCkteWs) — J3(h)A(hhhh) 



+ 



'1 / . 



{fhh 9khh + fhk 9l :i l 5 l 2 + fhh 9l 3 kh + fhh 9hkh + fhk 9h 



hh 



+ 



+ 



f hh 9hhkJ fhh 9mkhh fhk 9l 3 mhh fhh 9l 3 l 4 ml 2 fhh 9hmkh 

+fhk 9hhmh + fhh 9hhhm] (A^BtoW - B^A^) 

x m {fhh %khh + fhh 9hkhh + fhk 9l 2 hhk + fhh 9hhkh + fhh 9l 2 hkh 

X {A^ ^('a^kJsie) _ J^{h) A(hhkhh)^j 

x ™ {fhh 9khhh + fhk 9l 3 l s l 6 l 2 + fhh 9l 3 kkh + fhk 9l s khh 

+fhh 9l ikhh + fhk 9h hhh fhh 9hhkh fhh 9hhkh ) 

X f^d'a) gihkhh) _ J^(hh) A^hkhh) 



+ 



v m {fhk 9l 5 l 6 hh + fhh 9l 4 hhh + fhh 9khhh + fhk 9hhhh + fhh 9hhkh 
+fhh 9hkl 6 l 3 + fhk 9hhhh + fhh 9hhhk + fhh 9hhkh)_ {A^ 2 ^ B^ 1 '' 16 ^ t 



+ o( x 3 ) ■ 



(17) 
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B. Equations of Motion 

To determine the evolution to order C(x 3 ), we need the time derivatives of Xi(t), g%j(t), and 
gijkit). Take the commutator of products of generators with the Hamiltonian and subtract the 
disconnected parts to find: 



d 

Jt Xl ~ 



fij X a 



n ;j + [fij g ak + /;;, g ja 

+ [ x a (fij 9kl + fik 9jl + fil 9jk) + fij 9akl + fik 9jal + fil 9jka 
+ [ x a [fij 9klm + fik 9jlm + fil 9jkm + fi m 9jkl) 



Him 



~^~fij 9aklm + fik 9jalm + f H 9j 



kam 



fim 9j 



kla 



+ \%a [fij 9klmn + fik 9jlmn + fu 9jkmn + f i m 9jkln + fi n 9jklrnj ** 

+ o( x 3 ) 

J t 9^=[f?k9a 3 +f^a] 

+ [ x a (fik 9jl + fil 9kj + fjk 9il + fjl 9ki) + fik 9ajl + fil 9kaj + fjk 9ial + fjl 9ki 
+ \_ x a (fik 9lmj + fil 9kmj + fim 9klj + fjk 9%lm + fjl 9ikm + fj m 9ikl) 



;is) 



H (H) 



^ fik 9almj + /j2 9kamj + /i m fi'fciaj + /jfc 9ialm + /j7 9ikam + /j m fi'jfcic 
+ [^a (/jfe 9lmnj + /j/ 9kmnj + /j m 5tanj + /j n 9klmj 



£j(klr 



~^~fjk 9ilmn fjl 9ikmn ~\~ /jm 9ikln ~\~ /j n 9iklm 



■^j(klmn) 



(19) 



d 

+ [^a (/ ? 7 fi'jfem + fi'jfem + /fcj + fi'/jfe + // m &ifc + fkm 9lij) 

aj 9km "T 'y'afc S'jm ) + /jK&a fi'fcm + 9im 9ak) + f'kl(9 ia 9jm 9im 9aj ) 



+ fim(9lj 9ak + 9lk 9aj) + fj m {.9li 9ak + 9la 9ki) + fkmi.9li 9ja + 9lj 9ia) 
+ [ x a (fu(9mj 9nk + 9mk 9nj) + fimWj 9nk + 9lk 9nj) + firS,9l'j 9mk + 9lk 9mj) 
+fjl(9im 9nk + 9in 9mk) + fj m (9il 9nk + 9in 9lk) + /"„ (<7iJ 9mk + 9im 9lk) 
~^~fkl \9im 9jn + 9in 9jm ) + fkm (9 il 9jn + 9 in 9jl) + fkn (&i 

9jm 9im 9jl 



(20) 
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Recall that, through third order, there is no difference between the subtracted and connected 
correlators. Only the disconnected parts of the fourth order correlators appearing in Eq's. ([18]) 
and flT9|) are needed, since gijki = g^gu + QikQji + guQjk + C*(x 3 )- ^ equations only accurate to 
0(x 2 ) are desired, then all terms in Eq's. (p!8|)-(p0|) involving third (or higher) order correlators, 
as well as products of second order correlators, may be dropped.[] 

Given these equations of motion for the connected expectations of generators, one can use 
(|l2l) to describe the dynamics of any classical operator in terms of its symbol. If V = V{{xi}) 
is a (time- independent) function of the generators, then its time-dependent expectation value, 
at next-to-next-to-leading order, is given by 

(V(t)) ={V + 9ij (t) V<& + g ijk (t) + g ijkl (t) _ + 0( x 3 ) , (22) 

where Xi(t), gij(t) = s?j(t), and gijk{t) = Sijk{t) are to be obtained by integrating Eq's. (|18|)- (|20|) 
forward in time, using g ijk i(t) = fi^-(t) g k i{t) + g ik {t) gji{t) + gu{t) g jk (t) + C(x 3 ). 

C. Error Accumulation 

To any given order in x, we have a system of non-linear, first-order, ordinary differential 
equations. Initial conditions are imposed by specifying Xi(t=0) = (u\xi\u) and s^(t=0) = 
(u\xj ■ ■ ■ xi\u) conn , with \u) some chosen coherent state. Since s W(t=0) is C(x fe_1 ), and the 
equations for s^ k \t) involve only terms of order x k ~ 1 an d higher, we still formally have sV°'(t) = 
0(x k ~ 1 ) for t > 0. However, as the truncated equations of motion are integrated forward in 
time, errors accumulate; it is important to understand the rate of growth of this truncation 
error. 



7 The resulting next-to-leading order equations are simply 

j t Xi= ■>;,) ffW) + (/• g ak + a 9ja) + ,,, (/« gkl + /a gjl + f a ;hi) H (JH) + 0(x 2 } ; (21a) 

and 

j f 9ii = (ftk 9aj + f% 9l a) + x a (/& 9jl + ft g kj + f% ga + ft g kl ) + 0{ X 2 ) . (21b) 
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We are dealing with a system of equations which we can write as yi = Fi(y) + Gi(y) where 
{yi(t)} are the variables in our problem (that is, the x^s and sf.^'s), F(y) represents the terms 
we keep, and G(y) stands for everything thrown away by the truncation. Let yo(t) be the 
solution to the above equation with G = 0, and solve perturbatively, y(t) = yo(t) + e(t) with e 
small. Linearizing about y{t) = yo(t), we have 

e = f{t)e + g(t), (23) 

where //(£) = dFi(y (t))/dyj, gi(t) = Gi(yo(t)), and we have dropped C(e 2 ) terms. This 
linearized system of equations is easy to solve (at least formally). For t > 0, 



e(t) 



Te/o^^i c(0)+ f \Teti f ^ dt " 
Jo . 



g{t')dt'. (24) 



Here, T denotes time ordering (with smaller times on the right). If f(t) and g(t) are globally 
bounded during the time evolution, ||/(i)|| < /, ||<?(0ll — 9i "where || • • • || is some appropriate 
norm, then a crude estimate of the deviation of the true solution from the approximation is 

l|c(*)ll<e^||e(0)|| + sf(e^-l)//. (25) 

Of course for t small, errors grow linearly and ||e(t)|| < Il e (0)|| (1+ft) + gt + 0(t 2 ); with a 
truncation good to order x fc at t = 0, both ||e(0)|| and g will be 0(x k+1 )- 

In a general treatment, it is hard to do better than the crude bound (p5|) . In dynamical 
systems with only a few degrees of freedom, there typically are "regular" portions of phase 
space where perturbations grow only linearly with time |17j . However, it is not at all clear that 



this is applicable to the truncated quantum dynamics represented by Eq. (p3[). 

In simple examples discussed in the following section, we will find that for times of order 
X -1 / 2 , the shape of the wavefunction of the evolving state becomes so distorted that the formal 
hierarchy of correlators, ~ 0{x k ~ 1 )i upon which the truncation scheme is based, completely 
breaks down. In terms of the underlying quantum dynamics, if one considers the projection 
of the initial coherent state wavepacket onto the exact eigenstates of the Hamiltonian, what 
is happening for sufficiently large time is that the contributions of different eigenstates have 
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decohered to such an extent that the wavepacket has spread beyond recognition. Except for 
special non-generic cases (such as the harmonic oscillator, where there is no dispersion) one 
should always expect such decoherence to eventually set in. 



V. EXAMPLES 

We will discuss two examples of theories to which the preceding general results may be 
applied: the usual semi-classical limit of point particle quantum mechanics, and the large N 
limit of O(N) invariant vector models. For brevity of presentation, we will display explicitly 
only the first corrections to the leading classical approximation, but we emphasize that it is 
completely straightforward to include yet higher order corrections, such as the 0(x 2 ) terms 



displayed in Eq's. flT|)-(|20|). 



A. h — ► Quantum Mechanics 

Consider ordinary point particle quantum mechanics, in one dimension for simplicity. The 
coherence group G is the Heisenberg group, generated by {e{\ = {x/h,p/h, 1/h}. The formal 
parameter that controls how close the theory is to the classical limit is, of course, x — h- The 
rescaled generators of the coherence group, X{ = hei, include the position x and momentum p 
operators whose expectations will serve as classical phase space coordinates. The Heisenberg 
group, acting on a fixed Gaussian base state, generates conventional coherent states {\p, q)}, 
with wave functions given (up to an overall phase) by 

(x\p, q) = {-nh)^ 1 ^ exp j — ipx — \[x—q) 2 j. (26) 

We have arbitrarily chosen units such that our Gaussian base state has equal variance in x and 
p. Consider a Hamiltonian of the typical form H = ^p 2 + V(x), where, for simplicity, we have 
set the particle mass to unity. The equations of motion are, of course, 

-x=p, J t P = - V ^- (27) 
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+ 0(h 2 ) 
+ 0(h 2 ) 
+ 0(h 2 ) 

9pp 



We are interested in the time evolution of x(t), p(t), and the connected correlators g xx {t), 
9x P {t) = 9p X (t)i an d 9pp{t)i a ^ to order h. From equations (|1"8| ) and ( |I9"D we find: 

x = p + 0(/i 2 ) 

p = -V - + o(h 2 ) 

9xx 9xp 9px 
\9px) 9xp 9pp ^ 9 xx 

9pp ^ {9xp 4" fi'px) 

subject to the initial conditions x(0) = x , p(0) = po, ^ 2 
that to this order, det = g xx 9pp~ 9xp9px — 0(h 3 ) is a constant of the motion, and Eq's. (|28|) 
are equivalent to a Gaussian variational ansatz jjTj| (where one approximates the wave packet 
by a Gaussian with a time-dependent centroid and width). However, if we went to next-to- 
next-to-leading order in h it would become clear that our setup is different. For positive times, 
higher moments will not be given by simple algebraic expressions in terms of the centroid and 
variance, and the details of evolution will depend on the shape of the potential.^ 

As a trivial warm-up, consider the harmonic oscillator of unit mass and natural frequency 
Q: H = \p 2 + |fi 2 x 2 . The solutions to fl28|) are 



\% 9x P = ~9px = Notice 



x(t) = xo cos fit + (po/fi) sin fit , 
p(t) = po cos fit — (xq fi) sin fit , 



9xx(t) 

gpxtt) = 9xp(t) 

9 P p(t) 



cos 2 fit + fi -2 sin 2 fit 
i + (fi -1 — fi) cos fit sin fit 
cos 2 fit + fi 2 sin 2 fit 



(29a) 
(29b) 
(29c) 
(29d) 
(29e) 



Because the potential is quadratic these are exact. Equally simple is an inverted harmonic 



8 In our Gaussian initial state, all connected correlators higher than second order vanish at time zero, 
s( fc>2 )(0) = o. But these moments cannot remain zero unless the potential is harmonic. For example, us- 
ing Eq. ( pp| ) we find that s xpx = s ppx + s xpp — V" s xxx — V" (s xx ) 2 +0(h 3 ), showing explicitly that any nonzero 
V" will drive the skewness moments Sj,fc(t) away from zero. 
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oscillator. If one takes the Hamiltonian to be H = ^p 2 — |f2 2 x 2 , then the solution of the 



moment equations (p8 ) becomes 



x(t) = x Q cosh fit + (po/fi) sinh fit , 
p{t) = po cosh fit + (xq fl) sinh fit , 



9xx(t) 



cosh^ fit + fT 2 sinh z fit 



g * x (t) = g xp (t) = | \i + (fl 1 +fl) cosh fit sinh fit 



9pp(t) 



cosh 2 fit + fl 2 sinh 2 fit 



(30a) 
(30b) 
(30c) 
(30d) 
(30e) 



In both of these examples, the time-evolution of the variances are independent of xq and po. 
As one would expect, they oscillate (with twice the natural frequency) in the case of the simple 
harmonic oscillator, and grow (exponentially) for the inverted oscillator. 

As a more complicated example, consider the problem of small oscillations in a weakly 
anharmonic potential^] V[x) = \x 2 + [3x A . The moment equations (2^) become 



x 



(31) 



-x - 4/?x 3 - 12/? x g xx +0{h 2 ), 

9 xx = 9x P + 9 P x +0(h 2 ) , 

(g P x)* = 9x P = 9pp - 9xx - 12/5 x 2 g xx + 0(h 2 ) , 

g pp = -(1 + 12Px 2 ) (g xp + g px ) + 0(h 2 ) . 

We will solve these perturbatively; the two small parameters are j3q 2 and (5%. We will work 
to first order in (3fi (since we have omitted 0(h 2 ) terms in the moment equations), and will 
display explicit results through second order in f3q 2 . In principle, one could work to any order 
in /3q 2 desired. 

In order to keep our error estimates simple, we will treat the time as 0(1) (in units where 
the natural frequency is unity). This means we need not worry about the appearance of secular 



9 We choose the curvature of the potential at the minimum to equal unity, so that our chosen coherent 
states have the natural width for the unperturbed potential. This ensures that the resulting dynamics (such as 
oscillations of g^) are not merely reflecting purely harmonic oscillations. 
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terms — terms which grow as powers of t — and may solve Eq's. (|3T|) strictly perturbatively in 
the naive fashion. A straightforward calculation, with the initial conditions x(0) = q, p(0) = 0, 
9xx = Qpp = an d 9x P = 9p X = jih, leads to the solution 



x(t) = q jcost + (§/V) 



+ (ph) 
+ ((3k) (§/V) 



cos 3t — cos t — 12t sin t 

cos 5t — 24 cos 3t + 23 cos t + 96 1 sin t — 36 1 sin 3t — 72 1 2 cos t 
3 1 sin t 



f (cos 3t - cos t) - 18 1 sin 3t + 93 1 sin t-54t 2 cos t 



+ (/5g 2 ) 3 +0 (Ph)(f3q 2 ) 2 +0 (Ph) 



(32) 



with 



9xx{t) 

9x P {t) 
9 P p{t) 



\Ti{l-3f3q 2 
\%{\ + 3[3q 2 



1 — cos It + t sin 2t 
3 sin 2t + 2t cos 2t 
1 — cos 2t + t sin 2i 



+ o[(p q 2 ) 2 ]+om}, 

+ 0[((3q 2 ) 2 } + 0((3h)} , 
+ 0[((3q 2 ) 2 ] +0(J3h)}. 



(33) 
(34) 
(35) 



Examining the secular terms in Eq's. (|32j) and (|33|) , one sees that terms of order [3 k are accom- 
panied by at most k powers of t. This is a general result. It implies that our stated condition 
that the time be 0(1) is needlessly restrictive. For small (3q 2 and /3h, the perturbative expan- 
sions (|32"D and ([33|) are actually valid in the wider domain \(3q 2 t\ 1 and \(3Tit\ 1, provided 
a factor of t is included with each factor of j3q 2 or (3h in the error estimates. 

It is instructive to compare this treatment with the result of a perturbative quantum me- 
chanical calculation. Using the brute-force approach of first finding perturbed eigenstates and 
energy levels, and then evaluating the time-dependent expectation value x(t) by projecting the 
initial coherent state onto individual eigenstates and summing the resultant contributions, a 
rather tedious calculation using both wavefunctions and energies correct to O(0) leads to 

-,2 



x(t) = q e 



-(q 2 /h) sin 2 (3/37U/2) 



COS 



+ {/3ti) e-(<? 2 /^)sin 2 (3^ t /2)| cos 



3(3%) t + |-r sin(3/3M) 
q 2 

'l+3f3h)t + ^-sm(3(3ht) 
2n 



3 q 2 1 q 4 
2 ~ T ~ 32 ¥ 
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+ COS 
+ COS 
+ COS 
+ COS 



+ (Ph) e 



-(q 2 /h) sm 2 (9(3ht/2) 



COS 



2n 
q 2 

(l+9f3H)t+ ^sin(3^t) 
(1+12/3^) t + |rsin(3/3M) 
(1+15/m) t + |r sin(3/3M) 

3 (1+6/3^) t + ^- sin(9/3M) 



3 3 q 2 l? 41 

3 <t 
2 ft 
1 g 2 1 g 4 

1 g 4 ~ * 

32 n 2 
8 h 



(36) 



This result has 0(/3 2 ) errors due to the neglect of second (and higher) order corrections in both 
the eigenstates and energy eigenvalues. 

If one restricts t to be small compared to both l/|/3ft| and l/\f3q 2 \, then one may expand 
the result ( |36"D in powers of j3. Moreover, in this domain one may easily add in the leading 
secular 0[(j3h) 2 t] terms omitted from (p6|), which come from including the 0(/3 2 ) perturbation 
to energy levels while using unperturbed wavefunctions.^] One finds 



x(t) = gjcost + (§/3g 2 ) 



cos3t— cost— 12t sint 



+ (|/3g 2 ) 2 [96tsint-36tsin3t-72t 2 cost 



+ o 



—3 1 sin t 



{{3q 2 +{3hf 



+ (jm) 



-18 1 sin 3t + 93tsint - 54 r cost 



+ o 



{pq 2 +ph) 3 t 3 



(37) 



This result is perfectly consistent with the previous moment-hierarchy result (0), as it must 
be, except for the non-secular 0{(5 2 ) terms which are hiding in the first 0[(f3q 2 +f3%) 2 ] error 
term of fl37l) . If one includes second order perturbations to the eigenstates then these terms 
also coincide. 

In the semi-classical regime, where j3h <C /3q 2 , it is interesting to examine expression fl36|) 
when /3ht <C 1, making no assumption about the size of (3q 2 t. In this regime, the first, leading 
term of Eq. fl3"B"|) becomes 



10 This addition is q (t sint) [f± (/3q 2 ) 2 + ^ (/3q 2 )(/3h) + 18(/3?i) 2 ] . If one does not assume that fiht is small 
compared to 1, then including the 0(f3 2 ) energy shift in matrix elements of time-evolution operators unfortu- 
nately leads to an analytically intractable infinite sum for x(t). 
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x (t) = q e -!(^ 2 )W * 2 cos [(1 + 3f3h + \f3q 2 ) t ] + 0{q(3h) . (38) 

In other words, x(t) shows damped harmonic behavior, with a shifted g-dependent frequency, 
and with an amplitude which decays significantly on the time scale 



((3q 2 )(Ph)}~ 1/2 . (39) 



This implies that on this time scale, the initially well localized wavepacket has dispersed so 
much that its probability distribution is spread out over most of the classically allowed region.^ 
Hence, t^ should be regarded as a "derealization" or "decoherence" time. The higher order 
terms in Eq. (|36| ) all exhibit essentially the same behavior in this regime; each term oscillates 
with a (slightly different) frequency and has an amplitude which decays on the decoherence 
time scale td- 

Although it will have no bearing on our discussion, it is interesting to note that on yet longer 
time scales, when t is near 2n/{?>(3h) or integer multiplies thereof, the exponential factors in 
Eq. (|36| ) return to near unity, implying that the time-dependent state has "reassembled" itself 
into a recognizable wavepacket oscillating in the potential.^ Presumably, this is a reflection of 
the fact that this is an integrable single degree of freedom system. 

The existence of the decoherence time scale (|39|) has important consequences for the utility 
of any truncated moment expansion, such as Eq's. (|18|-f20D. If the wavepacket has spread to such 



n Of course, the fact that the amplitude of oscillations in the mean position x(t) decays on the decoherence 
time scale td cannot mean that the wavepacket has come to rest at the bottom of the potential while remaining 
a well-localized wavepacket, as this would violate energy conservation. In the semi-classical regime under 
discussion, the position of the initial wavepacket is significantly displaced from the minimum of the potential, 
q 2 ^> Ti, implying that the total energy is large compared to the zero-point energy. Therefore, a negligible 
mean position at large times necessarily indicates that the wavepacket has spread so much that its probability 
density, at any late time, is delocalized over the entire classically allowed region, and no longer "sloshes" back- 
and-forth to any significant extent. Within the classically allowed region, energy conservation implies that the 
wavefunction must have substantial variations on scales far smaller than the (square root of the) variance in 
position — which will be comparable to the width of the classically allowed region. 

12 Whether this "reassembly" persists in the exact solution, or is an artifact of our first order perturbative 
result, is not entirely clear to us. 
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an extent that it is significantly sampling all of its classically allowed region, while necessarily 
retaining structure on smaller scales, then the formal hierarchy of connected correlators, ~ 
fr*"' -1 , will have broken down. Higher order moments will not be small compared to lower order 
ones. Consequently, the moment expansion presented in the previous section can only be useful 
for times which are small compared to the decoherence time t^. 

The 1 / Vh dependence of the decoherence time (|39|) may also be seen in another very simple 
example. Consider the free evolution of a coherent state in the absence of any potential. As 
is well known, the width of the wavepacket grows without bound. The evolution equations 
for x and p are, of course, trivial, p{t) = p(0), and x(t) = x(0) + p(0) t. Hence, x 2 {t) = 
x 2 (0) + [x(0)p(0) + p(0)x(0)} t + p(0) 2 t 2 , and so for our initial Gaussian coherent state (with 
equal variance in x and p) , 

9xx(t) = g X x(0) + [g xp (0) + a»(0)] t + g PP (0) t 2 

= \h(l + t 2 ). (40) 

Here also, we see that for times of order hT 1 ^ 2 the hierarchy of correlators s^ k \t) ~ 0{Ti k ) no 
longer holds. 

We believe this to be a general result. Whenever a semiclassical system exhibits dispersion, 

-i in 

the decoherence time is expected to scale as n , and truncations of the moment hierarchy 
equations (|18H20|) will only be accurate for times small compared to the decoherence time. 



B. Vector Models 

Consider an O(N) invariant theory whose fundamental degrees of freedom form O(N) 
vectors. For simplicity, we will assume that the degrees of freedom are all bosonicQ and divided 
into a set of canonical coordinates {x l a } and corresponding canonical momenta {p^}. Here i, j = 
1,---,N are O(N) vector indices, while a,/3 = l,---,m distinguish different O(N) vectors. 



13 Extending this discussion to U(N) invariant fermionic models is completely straightforward. 
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These basic operators are assumed to satisfy canonical commutation relations, normalized 
such that [x l a ,fP^\ = (i/N) 5^ 5 a/ 3. In other words, we have chosen to scale both coordinates 
and moments by 1/yN compared to their textbook form. The small parameter controlling 
the approach to the classical limit is x = V-^i ^ has been set to unity. The Hamiltonian is 
assumed to be O(N) invariant, and we will completely restrict attention to the O(N) invariant 
sector of the theory. Consequently, the relevant Hilbert space Hn is the space of all O(N) 
invariant states, and all physical operators can be constructed from the basic bilinears 

N 

4=1 
N 

Ba P =T H*«ft+PK}. ( 41b ) 
4=1 

N 

&0 ( 41C ) 
4=1 

It will be convenient to regard x l a and p l a as the components of m x N matrices, so that the 
basic bilinears (O) may be assembled into m x m matrices, 

A = xx T , B = xp T - | i , and C = pp T . (42) 

Viewed as matrices, A and C are symmetric, while B is non-symmetric. The individual com- 
ponents of A, B, and C are all Hermitian operators acting on TIn- 
We will take the Hamiltonian to have the general form 

NH N = N[ltTC + V(A)]. (43) 

The overall factor of iV (given our scaling of coordinates and momenta by 1/yN) is exactly 
what is needed to ensure that the N —>■ oo limit is a classical limit in the framework of section II. 
The potential energy function V(A) may be any chosen scalar- valued function of a symmetric 
matrix A. The kinetic energy takes the simple form |tr(7 = \j2i,a(P a ) 2 ^ an degrees of 
freedom are scaled to have unit mass. Two specific examples in this class of models are: 

i) A single particle moving in a central potential in iV-dimensions. This is the simplest possible 
example; the theory has only a single O(N) coordinate vector [i.e., m = 1]. The Hamiltonian 
is 
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NH N = N 



~p- p + V(x ■ x) 



N 



\c + v(A) 



(44) 



where V{r 2 ) is now a function of just a one variable. 
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ii) An 0(iV)-invariant 4 field theory. The theory, defined on a spatial lattice, has field opera- 
tors (f) l s and conjugate momenta 7r*, where s labels the sites of some rf-dimensional lattice. The 



canonical commutation relations (after scaling <p and ft by 1/yN) are 
and the quantum Hamiltonian is 



{i/N) 5» 5 S 



NH 



~ 4 



A = • \ ; _ - k s ■ K s t- \ V0 S • V0 S + i/i 2 (f) s ■ <i) s + j (<j) s ■ 4> s ) (45a) 

(45b) 



{\ c ss + 1 [(-v 2 + ^)A SS ] I + 1 (i ss ) 



[Here V is a lattice forward difference operator, dot products denote the implicit sum over 
O(N) indices, and factors of lattice spacing are suppressed for simplicity] The number m 
of O(N) vectors [or the dimension of the matrices A, B, and C\ equals the total number of 
lattice sites. Ignoring the obvious notational changes (x — > 0, p — > tt), this theory has precisely 
the stated form of Eq's. (H^) . The lattice theory may, of course, be viewed as a natural 
discretization of the formal continuum theory where the field operators 4> l (x) and n l (x) depend 
on continuous spatial coordinates and 



NH N = N / (d d x) 



\ tt(x) ■ fc{x) + \ V0(s) ■ V0(x) + W }{x) ■ }(x) + \ U(x) ■ <j>(x) 



(46) 



Returning to the general discussion, a straightforward calculation shows that the commu- 
tators of the basic bilinears are 



— [A a /3, Arys) — -j [C a /3, C-ys] — , 
7" [Aa/3, B.ys] = A al 5/3$ + As 7 $a8 , 



(47a) 
(47b) 



In terms of coordinates and momenta which have not been rescaled by N 1 / 2 , NHn — jp 2 + N V(^x 2 ) 
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7" [B a /3, B jS ] = B lj3 5 a s — B aS 5 7/3 , (47c) 

Y [A a /3, C 7 s] = B ai 5/3S + _Bg 7 S a s + B a s <5g 7 + Bps S ai , (47d) 

Y [B a p, C-ys] = C/3 7 5 a s + Cps fia-y • (47e) 

In other words, the commutators of A, B, and C (as well as just A and B) close and these opera- 
tors generate a Lie algebra.0 The appropriate coherence group which will create suitable O(N) 
invariant coherent states may be taken to be the group generated by ( ant i- Hermit ian linear 
combinations of) the operators and {B a p}. Enlarging the coherence group by includ- 

ing the C a p operators among the generators is equally acceptable, but unnecessary. The group 
generated by {A1/3} and {B a p} alone satisfies all the conditions for producing an over-complete 
set of coherent states which behave classically as N — > oo. Including the C a p operators among 
the generators enlarges the coherence group, but has no effect whatsoever on the resulting 
manifold of coherent states. 

Acting on an initial Gaussian base state, the coherence group generates a set of coher- 
ent states where z is a complex symmetric m x m matrix, with positive definite real 
part, which may be used to uniquely label an individual coherent state. The position space 
wavefunctions of these coherent states are given by 



ifrJx) = det 



N / 



AT/4 



exp (-\NtYx T zx) . (4£ 



It will be convenient to decompose the matrix z into its real and imaginary parts by writing 

z = tjGT 1 - icu , (49) 
so that a = (z+z*) -1 and to = ~(z—z*). Both a and to are real symmetric matrices, and a 



15 The Lie algebra structure constants are s 



5rrS + &up&ij,6$a<y)) plus those trivially related by antisymmetry; all others 



vanish. The resulting Lie algebra of operators {A(a, b, c 
to the sp(2m) algebra represented by the 2m-dimensional matrices {A(a, b, c) 
etc., and a and c are symmetric. 



6„p6„ s 6 aj , and f Baf>CyS 



iN Y, a p( a ai3 Ap a +b a p Bp a +c a p Cp a )} is isomorphic 



}, where b = \\b. 



■ap\ 
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is positive definite. Using the fact that p l g\z) = ix l a z a/ 3\z), a short exercise shows that the 
coherent state expectation values of the basic bilinears are 



i „-i 



A(z) = a , B(z) — cuv , and C(z) = lucllu + \ a 



(50) 



The variances of these operators in the coherent state \z) are0 



AA 
CC 

( n BA \* _ AB 

( n CA \* _ AC 
[y-yS^P J — ya/3,78 

( n CB Y — n BC 

\y~y6,a/3 J — y a /3,jS 



faery &/3S da8 ^87] ) 

-a^+ujcujj)^ + (5 + iuoa)^ (~ - 



j_ 

N 



Clary ( 4 ^ 



3 )/37 (| _ ^)a<5 



A? 



iV 



N 



a a~/ (5 — i auJ )/36 + G/3 7 — 

~ - iau) ai (~ - iau)p S + (| 

(|a" 1 +a;aa;) j g 7 



v 2 io,ui') a ^ ( ^ 



■a 1 +a;acij) j g ( 5 



+ {\- iau) aS ( 



(51a) 
(51b) 
+u>au>)/3 7 , (51c) 
(51d) 
(51e) 
(51f) 



Given our choice of Hamiltonian (53), the operator equations of motion for the basic bilin- 



ears are 



-jI A a p — B a /3 + Bp a , 



B 



a/3 



d_ 

dt 

d c 



C a (3 ~ 2Aq, 7 ~V r. 



\/3 ^^«7 Vryfj , 

/ 



-IB^a Vi R — 2V aj B^p 



(52a) 
(52b) 
(52c) 



Here, V' is shorthand for the variation of V(A) with respect to the symmetric matrix A, 

6V(A) _ , 



V' = 

a/3 c a — 2 



5V(A) 8V(A) 



5A 



a/3 



5A 



/3a 



(53) 



and is defined so that 5V(A) = tr (V'5A).^\ 



16 For aesthetic reasons, we set g^ys = 9A all B^ s = (A a ffB jS ) - {A al3 ){B lS ) 1 etc. 



17 Note that with this definition, the matrix variation V reduces to an ordinary variational derivative in the 
case of a single vector (m=l). 
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Applying the general results (|l"8D and (|i9|) [actually, only (|21|) is needed] to the case at 
hand, one finds in a straightforward fashion the following equations, valid to next-to-leading 
order in 1/N, for the time evolution of the expectation values and variances of basic bilinears, 

A a p = B a/3 + Bp a , (54a) 



B a p = C af s - 2A arj VL - 2g AA V"g - g^a A *v K'kwCt + ®( N 2 ) > ( 54b ) 



d 

n — _9 r \r> _ 9 r \r< — I n BA _i_ n AB \ \r" _ ( n BA _i_ n AB \ \r" 

jT. — ^ D T]cx V V [3 ^-DriP V V a \9ria,^u "vw,Va J V vP,V \9riP,iJ.v ^ ^^VP ) V 



- fft&e v i#»a + B vP + °( iV " 2 ) > ( 54c ) 

together with 

9aP,y8 = 9a(3^5 + 9(3^8 + 9aP,yS + 9aP,8y + 0(N 2 ) , (55a) 

jL n BB - n CB + « BC> - 9« AB I/' - 9n BA V 

^ 9a/3,~/5 — 9 a p,y8 ^ 9a/3,-y8 ^9ari,y& VP "ctP,yn v rjS 



_ ( n AB _|_ n BA \ A V" 

- (gficp + g B .Lu) A ir y^u + o(N~ 2 ) , (55b) 

— n cc - -2n BC V - ?n CB V - ?n BC V - ?n CB V 

fl£ 9aP,jS ~ *9r)a,y8 %/3 9aP,V7 v rj8 *9r)P,y$ v ar) *9aP,r]8 v yr t 

- (gi% + sffij) (B v p V^ v + B m V^ v ) + 0(N~ 2 ) , (55c) 

9aP,yS = 9ap,yS + 9pa,y8 + 9aP,y8 ~ ^9aP,jn K;5 — ^9aP,/j,i/ A yv ^8,^ + 2 ) ) (55d) 

AC _ BC , BC _ 9 AB yi _ 9 AB yl 

9ap,y8 ~ 9ap,y8 * 9/3ari8 z i/Q/3,r;5 t/ 7»? z i/a/3,r)7 

- 2fl^ (5^ VJ,^ + Vfl^) + 0(N~ 2 ) , (55e) 

^. „BC _ CC _ 2 AC yi _ 2 BB y, _ ? BB yl 

^ 9aP,-y8 ~ 9a/3,-y8 *9oa),y8 v t)P ^9af3,r]8 v yr) & 9 a p,7]y v r)8 

_ ( n A C i f fA \ a yi, 

\9/j,u,y8 ' 9y8,fiu J ^-otr) v r]/3,nv 

- (9i B ap + 9aL„) {B n 8 V^ v + B m V^ v ) + 0(N- 2 ) . (55f) 

Here V^ )jS = SA ^^ , etc. Of course, g B s A a/3 = (gif a s) and so on ; since the basic bilinears 
Aap, B a p, and C a p are all Hermitian. 

As they stand, the (truncated) moment equations ( |51| ) and (^) are highly redundant. 
This is because the operators A a/3 , B a/3 , and C a p are not independent when acting on the 
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0(N) invariant Hilbert space T~Cn- For many purposes, it is preferable to reduce the evolution 
equations to a smaller set of independent observables. To see the redundancy, it is convenient 
first to note that the actions of B and A on any coherent state \z) are related, 

(B + il) \z) = xp T \z) = xx T (iz)\z) = A (iz)\z) . (56) 

Hence, the coherent state expectation value of A~ l B is directly related to that of 

^\{A- 1 B)jz) = iz aP -l(z\{A- 1 ) ap \z). (57) 

In a similar fashion, the coherent state expectation value of C may be expressed as 

(z\C\z) = (z\pp T \z) = (z\(iz)* (xx T ) (iz)\z) = (z\(px T ) (xx T )~ 1 (xp T )\z) 

= {z\{B + ^ A- 1 {B + ^i)\z) . (58) 

As noted earlier in section II, quantum operators are completely determined by their diag- 
onal expectation values in the over-complete coherent basis. Consequently, the coherent state 
relations floTD and ( |5*8"D suffice to infer underlying operator identities. The left-hand side of 
relation ( |57| ) is not manifestly symmetric under interchange of a and /?, but the right-hand 
side is symmetric under this interchange. Because Eq. ( p7| ) holds for all coherent states {|^)}, 
if one defines 

then (|57| ) implies that Vt a p = VLp a , so that Q = \\Q a p\\ is a symmetric matrix. Moreover, 
using the the commutation relations (|47|), one may verify that Q a/ 3 is Hermitian. [Demanding 
Hermiticity is what determines the coefficient of the second term in (|59|).1 Similarly, relation 
(|58|) implies the operator identity 

C = (B + A' 1 (B + \1) , (60) 



The following discussion assumes that coherent state matrix elements of A 1 exist, which requires N > m+l. 
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showing that the operators {C a p} are not independent of A and B [when acting on O(N) 
invariant states]. Inverting the definition ([59]) to express B in terms of Cl, 



B = AQ-i(^)i, 



(61) 



and using this, plus the Hermiticity of f2, allows one to rewrite expression (|60|) for C as 



c = CiaCi + \{i-^) 2 a- 1 . 



(62) 



Hence, within the O(N) invariant Hilbert space, instead of working with the basic bilinears A, 
B, and C [totaling m(2m+l) distinct operators], it is sufficient to use only A and Cl [totaling 
m(m+l) distinct operators]. These operators are, in fact, canonically conjugate "coordinates" 
and "momenta" . A short exercise shows that 



A-ap, A^§ 







%N 



Qaj3, AjS — Setups + Sasbp"! ■ 



(63a) 
(63b) 



If the complex symmetric matrix z parameterizing coherent states is separated into real and 
imaginary parts by writing z = \oT x — iu [as in Eq. fl49|)l, then the coherent state expectations 
of the canonical operators A and Cl are just a and u, respectively, 



(z\A\z) = a, (z\Q\z) = us . 



(64) 



[The first equality was previously noted in Eq. 



Re-expressing the quantum equations of motion (p2|) in terms of the independent canonically 
conjugate operators gives 

d 



A = AQ + QA , 
at 



(65a) 
(65b) 



where the "effective" radial potential 



V eS (A) = V(A) + ^ftrA- 1 



(66) 
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equals the original potential energy augmented by a "centrifugal potential". 

One may directly evaluate the evolution equations for expectations and variances of the 
canonically conjugate operators A and SI, or equivalently (and rather tediously) rewrite the 
previous equations (|5^ ) and ([55]) in terms of A and SI. Either way, one finds 

j t A aP = (An + siA) aP + gtSnfi + 9«Lp + °( N ' 2 ) > ( 67a ) 

j f Sl a p = ~(Sl 2 + W«U ~ 9^p - gXa + 0(N~ 2 ) , (67b) 

together with 

^ n AA _ A A q i SIA a , SIA a , „A4 q 

+ 9# m Sl v s + 9ap,r,6 + 9t$ m A vs + 9ap,r,s ^ + °( N ~ 2 ) > ( 68a ) 
d f qa \* _ d An _ An „ oq . on 4 , aq 

- - <4! 7 „ - 2 (v*)^ ^ + o(n- 2 ) , (68b) 



9aPa5 ~ ~9[ 



^ i/a/3,75 — Haf3,ri6 ^hv 9aP,-yr) 1 V z I v eg)~/6,iu> 9a.p,p,v 



9%: iS Sl av - O - 2 <V + ^(iV- 2 ) . (68c) 



Initial conditions corresponding to a given coherent state |z) (with z = |a _1 — icu) are given 
by A(0) = a and Sl(0) = u, together with the variances 

9at,l8 9aP,^\ 1 / daSap^ + a al aps ffiaS&Pl + Sarffips] 

.9aj3,-y8 9aP,-yS / t=0 ^ V _ 2 [^'Vy + ^ory'W I [ a a5 a /37 + a a 7 a /35] 

The next-to-leading order evolution equations fl67] ) and ( |68|) are directly applicable to any 
bosonic O(N) invariant vector model, such as the 4 theory defined by (f45|), whose Hamiltonian 



C(A^- 2 ) . (69) 



has the general form ([43). The dynamics is encoded in as efficient a form as possible; one has 
dynamical equations for the m(m+l)/2 pairs of independent phase space coordinates (0), and 
their variances (j68j) . 

In the special case ( pH|) of a single vector (corresponding to a point particle moving in an 
N- dimensional spherically symmetric potential) one may drop all the indices and the next-to- 
leading order evolution equations become 
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jA = An + nA + g An + g nA + 0(N- 2 ), (70a) 

f f n = -n 2 - 2K' ff - gnn - 9aa Kff + 0(N~ 2 ) , (70b) 

^ 9aa = 4g AA n + 2(#m + ^ A + G(iV- 2 ) , (70c) 

^ GM* = J f 9An = 2gnn A-2g AA V£ + 0(N~ 2 ) , (70d) 

^ <?nn = -4#nn - 2(&u + <u n ) + 0(iV~ 2 ) , (70e) 



with initial conditions given by A(0) = a, Q(0) = u, and 

+ 0{N~ 2 ) . (71) 



9aa 9 au \ 2 ' " 2 



,9ua gnn J t=0 \ — 2 i a , 
From Eq's. ([70]) and ([71]) one may again see that to next-to-leading order, the determinant of the 
variance matrix on the left-hand side of ( [HP is a constant of the motion, detg^(t) = O (N~ 3 ). 
To this order, our method gives exactly same predictions as the Gaussian approximation of 
||18|| . One may, of course, systematically extend the treatment to higher order in 1/N simply 
by specializing the next-to-next-to-leading order results in section [TV| . 

The evolution equations QT0"| ) in this single-vector case may be cast in a more transparent 
form by defining radial position and momentum operators via 

A = r 2 , tt = \ipf- 1 +f- 1 p) , (72) 

or equivalently 

r = A 1 ' 2 , p = A l ' 2 Q-^A- 1 / 2 . (73) 
These operators are canonically conjugate, 

i\p,r] = l/N, (74) 



and a short exercise rewriting the quantum equations of motion (p5|) yields 

±r = p, (75a) 
±p=-U'«(r), (75b) 
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where 



U pff (r) = K ff (r 2 ) - 



8N 2 r 2 

= ^ 2 ) + |(l-f)(l-i)r- 2 , (76) 

and U' c q = dU c s/dr. This is a well-known result: s-wave dynamics in an iV-dimensional central 
potential is equivalent to one- dimensional quantum dynamics in an effective radial potential 
U c g containing an additional "centrifugal" potential ^g^i^T^ which is non-vanishing in all 
dimensions other than 1 and 3 [21,22]. As seen in the commutation relations (|74|) , the parameter 



1/N plays the role of h so that the large N limit is precisely equivalent to the semiclassical 
limit of ordinary one-dimensional quantum mechanics. 

The next-to-leading order evolution equations ( ]70D for the coherent state expectation values 
and variances of A and Q may be easily be converted to equivalent next-to-leading order 
equations for expectations and variances of p and r. One finds, 

±r = p + 0(N- 2 ) } (77a) 
j t P= -U' eS - \g rr U% + 0(N~ 2 ) , (77b) 

4 9rr = 9rp + 9pr + 0(N~ 2 ) , (77c) 
J t (gprT = J t 9r P = 9 PP - 9rr + 0(N~ 2 ) , (77d) 

j t 9 PP = ~(9r P + g pr ) U'J S + 0{N- 2 ) . (77e) 

Through next-to-leading order, these evolutions equations are identical to the evolution equa- 
tions (|28| ) for the usual semiclassical limit.0 The initial variances differ, however, due to the 
differing shapes of the initial wavepackets ( ^6]) and ([!§). For our large coherent states, 

r 2 pr+i \ 

+ 0(N~ 2 ) , (78) 

pr—i p +r / 



/ 9rr 


9r P 




\9pr 


9 PP l 


2N [ 

t=Q > 



19 This equivalence persists to all orders, of course, reflecting the exact correspondence between the operator 
equations of motion (|2^) and (f75|) . 
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[and once again det g( 2 \t) = O (iV~ 3 )]. The form of this variance matrix (including, for exam- 
ple, the growth in the variance g rr with increasing r) reflects the fact that the underlying O(N) 
invariant coherent state wavefunctions are not constant width one-dimensional Gaussians, but 
rather iV-dimensional Gaussians centered at the origin with variable width. Hence, the posi- 
tion of the peak in the resulting radial probability distribution is positively correlated with the 
width of the radial probability distribution about this peak. 

For any given choice of the potential, one may integrate the five equations ( |77| ) forward in 
time and obtain results which are accurate to 0(N~ 2 ) [for times of order unity]. For better 
accuracy, one could extend the treatment to include higher order correlations, as detailed in 



section [IV . 

In light of the above exact correspondence between the 0(iV)-invariant dynamics of the 
single- vector model fl44]) , and ordinary one-dimensional quantum dynamics in the the effective 
radial potential (|76|) with N playing the role of h, the previous discussion of stability of the 
truncated moment equations in the semiclassical limit immediately carries over to the large N 
dynamics of the single-vector model. In particular, this means that one should expect to see a 
decoherence time which scales as N 1 ^ 2 , beyond which truncations of the moment hierarchy are 
no longer useful. We have no reason to believe that the scaling of the decoherence time with N 
will be different in more general vector-like large N theories, such as the 4 field theory ([45]), as 
compared to the single-vector model. Although we have no compelling proof to offer, we expect 
that a decoherence time of order N 1 / 2 is a generic feature of vector-like large TV theories.F] 



20 It is interesting to note that, in contrast to the previous discussion of the semiclassical % — > limit, examining 
iV-dimensional free motion in the absence of any potential does not provide an example illustrating breakdown 
of the moment hierarchy based on O(N) invariant coherent states. This is because the growth in the width 
of a spherically-symmetric Gaussian wavepacket is perfectly represented by a single one of the variable-width 
O(N) invariant coherent states (|4|), unlike the earlier situation with fixed-width coherent states. Hence O(N) 
invariant free motion is highly non-generic. For O(N) invariant free motion (in the general case where z is an 
m x m matrix and Hn = ^trC), one may show that the exact time evolution maps an initial coherent state 
\zq) into another coherent state \z(t)) with z(£) _1 = z^ 1 + it 1. The operator equations of motion ( |52|) may also 
be integrated exactly and show that C(t) = C(0) is a constant of the motion, while B = B(0) + 2Ct, and A = 
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VI. CONCLUSIONS 



We have shown that a systematic hierarchy of time-local evolution equations for a minimal 
set of equal-time correlation functions may be derived in any theory having a classical (or 
large- N) limit which fits within the general framework of section ITJ. Truncating this hierarchy 
at the level of fc'th order moments (i.e., retaining up to /c-point connected correlators) yields 
results which are accurate up to order 1/N k . 

However, it is clear that the t — > oo and h — > (or N — > oo) limits are non-uniform. At least 
in simple one degree of freedom (or single vector) models, we have argued that integrating the 
truncated moment evolution equations forward in time yields results which, generically, cease 
to be a good approximation to the true quantum dynamics beyond a decoherence time which 
scales as hT 1 ^ 2 (or y/N). The ordering of connected correlators which underlies the truncation 
of the moment hierarchy is only valid for times small compared to the decoherence time. Going 
to higher orders in the truncation scheme will not, in general, yield results which remain valid 
for parametrically longer time intervals. 

We expect, but have not demonstrated, that this y/~N scaling of the decoherence time is a 
general feature of large N quantum dynamics. It would obviously be worthwhile to investigate 
this further, particularly in large N models with many vectors. In, for example, an O(N) 
invariant lattice 4 field theory, it would clearly be desirable to understand the dependence of 
the decoherence time on the energy of the initial state and the lattice volume. If the yN scaling 
of the decoherence time is generically true this would, for example, imply that one cannot use 
truncated hierarchies of large N evolution equations (at least of the form considered here) 
to study the non-equilibrium dynamics of thermalization or hydrodynamic transport, as the 
relevant time scales for these processes scale like N in the large N limit [TjJ. We hope that 
future work will shed light on these issues. 



A(0) + | B(Q) + B(0) T t + Ct 2 . This implies, for example, that for large time the variance g^iuv ~ an d 
so grows without bound. However, the mean value (A) grows quadratically with t, and hence rms fluctuations 
remain bounded and of order N~ x / 2 for all times. 
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